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Abstract. The effects of macro- and microlensing on the spatial and 
temporal characteristics of images of remote sources, observed through 
the inner regions of lensing galaxies are discussed. A particular attention 
was given to the case, when microlenses, - stars, star-like or planetary 
bodies, - are situated near the critical curves of macrolenses, - galaxies, 
stellar clusters , etc. The investigation is of interest for the gravitational 
lens (GL) systems, where the lensed images are observed close to the 
critical curve of a macrolens. Annular, arched or confluent images should 
be regarded as an indication to such a proximity. Numerical simulation 
allowed to determine structures of critical curves and caustics, formed 
by macro and microlenses, and to evaluate possible distortions, caused 
by microlenses for various locations with respect to the critical curve of a 
regular lens. The difference of our results from those obtained earlier with 
the standard (linearized) approach to describe the regular gravitational 
lens was shown to be the larger the closer to the critical curve the mi- 
crolenses are situated. In addition to spatial redistributions in the visible 
images, complicated deformation of their light curves occurs in gravita- 
tional lensing of variable in time and extended sources. The magnitudes 
of the temporal variations depend not only on the GL parameters (e.g., 
mass distribution), but on the parameters of the source as well, such as 
linear dimension of the emitting region, its location with respect to the 
critical curve, and the impulse duration. The analysis shows, that in this 
case GL acts as a filter, which passes slow temporal variations without 
degradations, and smoothes the rapid ones. 



1. Introduction 

Uncorrelated variations of brightness in the images of a distant source, splitted 
by the gravitational field of the lensing galaxy are indicative of microlensing 
events. The phenomenon is clearly seen, e.g., in Q2237+0305, which is known 
to consist of four lensed quasar images within a 1-arcsec radius circle around 
the center of the lensing galaxy. The uncorrelated brightness variations of these 
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images indicate that the effects of microlenses are mainly observed, but not the 
intrinsic quasar brightness variations. Because of significant difficulties of theo- 
retic analysis, microlensing phenomena are investigated mostly with numerical 
methods. A standard procedure for numerical simulation consists in the fol- 
lowing. The effect of all microlenses, randomly distributed inside the galaxy, is 
described in a twofold manner. The microlenses, which are the most close ones 
to the observed image of the source, are specified with their actual deflection 
angles. The number of such microlenses is not expected to be larger than sev- 
eral hundreds, (Paczynski 1986; Seitz & Schneider 1994; Blandford & Narayan 
1992). The rest of microlenses, which are far from the line of sight, as well as 
the smoothly distributed matter, are specified by introducing a deflection angle, 
averaged over the interstellar space. The further simplification of the problem is 
due to linearization of the mean deflection angle of the macrolens near a selected 
lensed image, (Seitz & Schneider 1994; Blandford & Narayan 1992; Schneider 
et al. 1992). The procedure of the local simplification works well for isolated 
source images, which are essentially far from the macrolens critical curve. It is 
of interest, however, to investigate also an overall interaction of macro and mi- 
crolenses for the lensing galaxy as a whole. Therewith, a possibility appears to 
follow the changes of spatial density of microcaustics, and the velocities of their 
displacement depending on their location with respect to the critical curve of 
the macrolens. In these cases, typical spatial and temporal scales of microfocus- 
ings substantially differ from those obtained with the linearized approach. The 
analysis is of a particular importance in studying microlensing by a star, which 
is in its turn surrounded by a large amount of planetary-mass objects. Such a 
complex GL will be further considered as a system of macro-milli-microlenses. 

In addition to spatial redistribution of brightness in the apparent images, 
complex deformations of their light curves also take place in gravitational lensing 
of sources, which are extended in space and variable in time. The amount of the 
deformation depends not only on the parameters of the GL, - mass distribution 
in the lens, - but on the parameters of the source as well, such as the emitting 
region dimension, its location with respect to the critical curve, and on duration 
of the source brightness spikes. The amount of light curve deformations turned 
out to be different for different lensed images. GL can be regarded as a multi- 
channel filter, which passes slow brightness variations and smoothes rapid ones. 

The aim of this work was to analyze theoretically and simulate numerically 
the structures of critical curves and caustics of complex gravitational lenses, 
which consist of macro- and microlenses, for various locations of microlenses 
with respect to the macrolens critical curve. Also, peculiarities of gravitational 
focusing of the source with the arbitrarily distributed masses were investigated 
for the case, when the source is extended in space and variable in time. To obtain 
numerical estimates, the results were applied to the Einstein Cross Q2237+0305. 

2. Spatial changes 

A complex GL can be presented as a sum of two constituents, with one of 
them being due to the masses, contained in compact objects, such as stars, 
star-like and planet-like bodies, and the second one, connected with a diffuse 
matter, such as dusty and dusty-gazeous clouds. Compact objects, referred to 
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as microlenses (or millilenses) hereafter, can be defined as point masses Mj. A 
total amount of microlenses n in the lensing galaxy is sufficiently large. For 
example, N « 10 11 -=- 10 12 for a spiral galaxy similar to Q2237+0305. The effect 
of microlensing was at first analyzed for a single microlens, and the results were 
then generalized to the case of many microlenses. Let us introduce a Cartesian 
coordinate system XYZ with the origin in the galaxy center of masses, and the Z 
axis crossing the point of observations, (Fig. 1). Let Z p and Z s be the distances 
between the GL, the observer P and the source S, respectively. Since Z p and 
Z s are much larger than the size of the region, where the light beams undergo 
deflection, we consider the GL as a thin phase screen, (Bliokh k, Minakov, 1989), 
which is superposed with the Z = plane. Having selected a single microlens 
from a totality of N microlenses, located at the lens plane and characterized 
by the impact parameter ip = if) m , we have for the deflection angle at the lens 
plane, (Z = 0): 

Here, tp = (V'xiVfy) is the angle of observation, tpi a coordinate and r g i = 
2GMi/c 2 a gravitational radius of the i-th microlens, G is the gravitational 
constant, and c is the velocity of light. ©^/(VO is a component vector of the 
deflection angle, connected with the diffusely distributed matter. Supposing a 
multidimensional distribution function of random quantities r g i,ipi and a distri- 
bution law for the diffuse matter in the galaxy to be known, one can obtain the 
average defection angle < > and the variance < (<50) 2 > from Eq. 1, (see 

e.g., Bliokh & Minakov 1989, Schneider et al. 1992). 

The GL with the average deflection angle < 0(VO > is called a regular 
lens, or a macrolens. The approximation suits a GL, which has the densities of 
stellar masses and diffuse matter " smeared" over sufficiently large regions of the 
interstellar space. 

Consider a "small" region of space, located near the m — th microlens, 
which has the size, from one hand, small enough as compared to the typical 
interstellar distance *$> s t in the particular region, and from the other hand, is 
larger than or comparable to the characteristic angular size of the Einstein ring 

= ^2r gm Z s /Z p (Z p + Z s ) for the given microlens, (* m < \if - if m \ < * st ). 

Accounting for the smallness of the RMS deflection angle (Bliokh & Minakov 
1989, Schneider et al. 1992), one can easily show, that in fact, the second and the 
third addends in Eg.(l) do not differ from the mean deflection angle < Q(tp) > ■ 
With all the above noted taken into account, the equation of the complex 
lens can be written as: 

& = + # < eW > = HA (2) 

where ip s is the actual location of the "point" source, and Z = Z p Z s /(Z p + Z s ). 

The lens equation can be better written in normalized quantities V = 
il)/^l,v s = ^ s /^i,v m = ipm/^L- It is convenient to determine the normalizing 
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angle $l as the angular radius of the Einstein ring, which would be observed 
for the given macrolens, provided it did not contain an asymmetric constituent 
of the mass distribution. In the normalized quantities, the equation of lens (2) 
can be presented as 

v s = F M (u) - Y m l~y = F(u), (3) 

where Fm{v) = ^+ < 0(z?) > is a constituent, connected with the macrolens 

galaxy, and T m = ^ 2 m /^f << 1 is a small parameter, that is approximately 
equal to the ratio of the microlens mass to the total mass of the galaxy. The 
gravitational lens system Q2237+0305 was selected as an example, for which 
various models of mass distribution are known, (see, e.g., Kent& Falco 1988, 
Minakov & Shalyapin 1991, Wambsganss & Webster 1997). Basing on the ob- 
servational data, let us choose the model to be simple enough, but in such a 
way, that it would be able to produce the observed cross-shaped structure. Let 
the regular constituent of the deflection angle < 0(z?) > be presented as a sum 
of two terms, with one of them, 0d(z7) being connected with a massive disc, 
and the other one, ® c {p), - with a compact nucleus. To simplify the analysis, 
suppose the nucleus and the disc are spherically and elliptically symmetric, re- 
spectively, with their centers coinciding. Superpose the directions of OX and 
OY axes with the axes of symmetry of the disc component. Since the typical 
scale for changing the surface density of the disc component in Q2237+0305 
exceeds noticeably the region, where four quasar components are observed, - 
about 2", - one may use the following approximation for @d(v) near the center 
of the galaxy: Q d (v) ps 9 d (0) + (z?V)0 d (O), where V = -£^e x + -£-e y , e x and e y 
are unit vectors of the Cartesian coordinate system. Having supposed charac- 
teristic size of the galaxy compact nucleus to be smaller than the region, where 
the four images are observed, let us approximate the spherical component by 
a point mass, equal to the mass of the galaxy nucleus Mq- With the accepted 
simplifications and suggestions, and having in mind, that for the chosen model 
of the disc 0^(0) = 0, the regular component of the macrolens equation can be 
brought to a form, (Minakov & Shalyapin 1991): 

Fm{v) = v{l - ^ ) + av x e y - av y e y . (4) 

Here, a is a dimensionless parameter that determines the amount of asym- 
metry, introduced by the disc component near the galaxy center. For the 
Q2237+0305 model, the normalizing angle ~ 0.9", T m can be estimated 
to be T m ~ 1CP 9 -T- 1CP 10 , and the asymmetry parameter can be selected equal 
to q « 0.16. Therefore, the inequality T << a << 1 holds. 

The equation of the critical curve of a GL V = V CT is obtained from the 
condition that Jacobian of transformation from variables V s to variables V is 
equal to zero: 

,-V) = ||^|=0, (5) 
where q(p) is the amplification factor. 
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If we substitute the solution of v C n 

found from (5), to the initial lens equa- 
tion (3), we will obtain the equation of caustic in the source plane: 

Vcs = F M {y C r) - r m _ " m ^ 2 - ( 6 ) 

\U cr Vrn) 

As the first step of the further analysis, let us determine the critical curve 
v = zT^f and caustic v cs = of the macrolens, having supposed T m = in Eq. 
(3). For example, for the model (4) in the polar coordinates V = (u, <p, ) in a 
linear approximation with respect to a small parameter a, we have: 

i^ / ( V ?)^l + -cos2 V 9, (7) 

with the azimuth angle counted off the X axis. 

The equation of caustic of the regular GL in the source plane in the same 
approximation is: 

u^ifp) = 2a cos 3 <p, vf sy ^) = -2a sin 3 <p. (8) 

The model of a regular GL under consideration is seen to produce an oval 
critical curve in the lens plane, with an approximately unit average radius, - 

« 1, while the corresponding caustic curve in the source plane is diamond- 
shaped, with the typical size of \v%\ « 2a. It can be easily shown, that four 
images will emerge in projection of the point source S inside the diamond caus- 
tics (8), with two images outside the critical curve, and two ones inside the 
curve (7) in the lens plane, (see, e.g. Minakov & Shalyapin 1991). Because of 
smallness of a, the images are located close to the critical curve (7). 

3. Effect of microlenses 

Regarding the equations of the macrolens critical curve and caustic to be known, 
one may investigate deformations of the curves for various locations of a mi- 
crolens v m with a specified parameter T m . A simple method was applied to nu- 
merically solve Eq.(5). The problem of finding the roots of equation q~ l {p) = 
can be substituted by a more simple one. Thus, it is easy to numerically de- 
termine the regions of positive, q~ l {v) > 0, and negative, q' 1 ^) < 0, values of 
q~ l (v). Excluding singularities, connected with the point mass approximation, 
the problem is reduced to determining the boundary between the positive and 
negative values of q~ l (v). The simulation has shown, that a behavior of criti- 
cal curves and the relevant caustics of a complex GL depends not only on the 
distance between a microlens and the critical curve of a regular GL, but also 
on the sector, where the microlens moves. Critical curves and caustics for the 
Q2237+0305 model were calculated for different locations of the microlens, - in- 
side and outside the critical curve of the regular GL, and exactly in it, - and are 
shown in Fig. 3a,b,c. The cases of the microlens displacement along the X and 
Y axes were examined. To ensure better view of macro- and microcurves behav- 
ior in the same picture, a sufficiently large value for V was selected, T = 0.001. 
The curves for T = 10~ 9 will be also presented below. 
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For a microlens outside the critical curve of a regular GL, - Fig.2a, q^ 1 (v) > 
region, - the curve v = v cr consists of two closed curves. One of them is a 
slightly deformed macrolens critical curve v^f , and the other one appears around 
the microlens. The similar behavior of the critical curve of a complex GL is ob- 
served also when the microlens is inside the critical curve of a regular macrolens, 
- Fig. 2b, < 0. In contrast to the case 2a, however, the microlens critical 

curve is not a single deformed oval around the microlens (Einstein "ring"), but 
consists of two closed curves on either side of the microlens. The most inter- 
esting case, when the microlens is situated right at the macrolens critical curve, 
is shown in Fig. 2c. The critical curves of macro- and microlenses merge into a 
joint continuous curve in this case. 

A numerical analysis allows to clearly demonstrate a complex deformation 
of critical and caustic curves, but mathematic analysis is needed to obtain qan- 
titative estimates of the expected effects, that would make it possible to follow 
their dependence on the model parameters. Determining of the real roots of 
Eq.(5) for the given values of displacement u m and parameter T m << 1 is the 
major task of the analytic investigation. 

The following conclusions have been made in the present theoretic analysis, 
applied to a particular case of Q2237+0305, (Minakov & Vakulik 2000): 

• A microlens distorts a critical curve and caustic of a regular GL. 
The amount of deformation is proportional to T m and decreases rapidly 
as the distance between the microlens and the macrolens critical curve is 
growing, see Fig. 3. Since rather large number of microlenses may be 
present near the critical curve simultaneously, their joint effect will cause 
"blurring" of boundaries of the critical curve and caustic of the macrolens. 
In Fig. 4, the structure of critical curves and caustics for Q2237+0305 is 
shown for the case of three hundred microlenses, located near the v^J. . 

• A structure and behavior of critical curves near the microlens depends not 
only on the distance to the critical curve of the regular GL, but on the 
sector of angles as well. Characteristic dimensions of critical curves and 
caustics are proportional to VT m . 

• The strongest deformation of critical curves and caustics, proportional to 

1/3 

Tm is observed when microlenses are situated very close to the critical 
curve of the regular GL. The curves for macro- and microlenses merge 
together to form complex continuous curves. 

• The difference between the results, obtained with the proposed approach 
and those of the linear approximation is the larger the closer a microlens 
to the critical curve, as can be seen from Fig. 5. 

• The effect of interacting microlenses increases noticeably in their approach- 
ing the critical curve of a regular GL. 

• If even microlenses are distributed homogeneously, spatial distribution of 
their microcaustics is substantially inhomogeneous near the macrocaustic, 
see Fig. 4. Such an inhomogeneity results in spatial scale of microlensing 
being dependent on the macroimage location with respect to the macrolens 
critical curve. 
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The theoretic analysis and the obtained estimates are capable of examining 
more complicated situations as well. Suppose, for example, that a complex GL 
system consists of a macrolens galaxy, millilens star, with the mass of order of 
Mq, and of a large number of microlenses, - planetary bodies with the masses 
of 10 -5 Mq. If a millilens is situated far from the critical curve of the macrolens 
galaxy, the effect of focusing can be more easily analyzed in the coordinate 
system superposed with the millilens. In doing so, the gravitational field of 
the galaxy may be considered in the linear approximation, while the effect of 
microlenses, randomly distributed around the millilens, should be considered in 
the way similar to that described above. It is only necessary to renormalize the 
angles to the quantity corresponding to the Einstein ring of the millilens. 
The parameter will be proportional to the ratio of the mass of the fi-th 
microlens to the mass of the millilens. 



4. Temporal changes 

A complicated redistribution of brightness in the observed image takes place in 
gravitational focusing. If parameters of the source and the lens do not change 
in time, the brightness distribution I p (tp), observed through the GL is related to 

the actual source brightness distribution I s (ip s ) by a convolution-like formula, 
(Bliokh & Minakov 1989): 



I s (ii s )6 - F(Y>) d^ s . (9) 

-oo 



Here, 5 ip is the delta-function. The formula represents mathematically the 

property inherent in all the nonabsorbing lenses to remain the brightness along 
a particular light ray unchanged. 

Several compact bright images (macroimages) , separated by a few arcsec- 
onds are usually observed in actual gravitational lenses. The macroimages have 
dimensions approximately coinsiding with those of the source emitting regions. 
For example, they are about 10~ 4 -j- 10~ 6 arcseconds for quasars in optical 
wavelegths. For radio wavelengths, angular dimensions of macroimages may 
be several orders of magnitude larger. Contemporary instruments, - telescopes, 
radio telescopes, - do not have a resolution high enough to examine the fine 
structure of macroimages I p (ip). Therefore, the source brightness, integrated 
over angles will be observed: 

/oo _^ 
i P W#. (io) 
-oo 

Amplification of the source brightness, resulted from gravitational focusing of a 
stationary source, is characterized by the amplification factor g, which is defined 
as: 

q = J p /4°\ (11) 

where Jp^ is the source brightness, which would be observed without the GL. If 
J p is ment to be a total brightness of the whole distorted image of the source, then 
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q determines the total amplification factor of the lens. If J p is the brightness 
of the i— the macroimage, then q means the amplification factor of the i— th 
macroimage, In this case, coefficients q and qi are associated by a simple 
relationship: 

i 

In the more general case, when the brightness of the source varies not only 
in angles, but also in time, - I s = I s (ip s ,t), - while the GL parameters are 
virtually unchanged, equations (9) and (10) can be generalized to take the form 
(Minakov & Nechaev 1998): 

POO p _,_-|r_ -> -i 

I p (if>,t) = I 9 il> 8 ,t-tp{il>,il> s ) 6 fa-Fty) dif>, (13) 

J — oo 

i P (4>,t)di>. (14) 

-oo 

tp(ip,ip s ) is a signal propagation time along a virtual ray, that connects an ar- 
bitrary point of the source (ip s ) with an arbitrary point at the lens plane (ip), 
and, further, with the observer P. Taking into account that the local velocity 
of the signal propagation in the weak gravitational field with a potential <I>, - 
|$|/c2 << 1, - is equal to c/n(f), and introducing the effective refraction index 

of the "medium" n(r), (Bliokh & Minakov 1989), the propagation time t p (ip,ip s ) 
can be determined by integrating n(f) along a curvilinear ray L, that associates 
arbitrary points of the source, the lens and the observer: 



tp 



1 J n(f)dl. (15) 



The effective refraction index can be representd as n(r) = 1 — 2$(r)/c 2 . 

Using a "thin" lens and paraxial optics approximation, intergration along 
the curvilinear ray in (15) can be replaced by integration along the lines, con- 
necting a point element at the source surface with a point element at the lens 
plane, and the lens point with the observer. After simple transformations, we 
have the following expression: 



rr i rr i Z p (Z p + Z s ) 1 1 Z p + Z s 2 1 ,2 77 \ 
Z p + Z s + 1 I- — i) s + -ip -i/npsj 



+tgr(lf). 



(16) 

The first addend in this formula is the propagation time for the signal in 
the empty space along the broken path " source- lens-observer" , see Fig. 1 . The 
second term, t gr , is due to the presence of gravitational field. For the weak 
gravitational field, t gr (ip) can be calculated by integrating of a gravitational 

addition 2<£/c 2 in (15) along the rectilinear rays. It is easy to show, that t gr {ip) 
can be reduced to a form: 

2G f AlZ 
V W = — / dipatf) In 8 + const, (17) 

c 6 Jw Z p (ip' - ip) 2 
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where cr(ip) is the given surface mass density distribution for the lens, [a] = 
kg/rad 2 , and integration is carried out over the whole surface of the GL in 
projection to the sky, and const is an insignificant addition, affecting only a zero 
point of the time. Eq. (16) and (17) hold for the low-redshift GL systems, when 
the cosmology factor can be neglected. 

Eq. (13) has a simple physical meaning. The point spread function of 

the lens, S tp s — F(ip) , associates every particular angle of observation with a 

certain "point" element of the source surface, having an actual angle ijj s = F(ip) 
and a surface brightness I s (ip s ,t). Temporal changes of I s (ip s ,t) will be recorded 
by the observer at an angle ip with a time delay of t p (ip) = t p (ip,^ s = F(ip)). 
The exact coincidence between the temporal variations of the source brightness 
and brightness of its images is seen to occur only for the point sources. In 
focusing extended sources, their actual light curves will differ from the light 
curves of the lensed images. 

To simplify the further analysis, we studied a variable source with the Gaus- 
sian brightness profile: 



Jo 



exp 



2*2 



f(t,4> s ). 



(18) 



where *o an d * s are angular dimension and angular coordinate of the source 
brightness maximum, respectively, and f{t,ip s ) is a dimensionless function of 
the order of unity, accounting for the temporal variations of brightness. It is 
supposed further, that the source varies in time as a whole according to: 



f(t^ s ) = ^-ex P l-^-Yh(t), 



(19) 



where T s is a characteristic duration of the brightness flash, h(t) = 1 at t > 0, 
and h(t) = at t < 0. 

This model refers to the source, which flashed at t = 0, and then is gradually 
fading, with a characteristic decay time T s . The Jo quantity in Eq. (18) is the 
full energy, emitted by the source: 



/oo rco 
dtps / dtl s (ip s ,t). 
-oo J — oo 



(20) 



For the given model (18), and according to (13), the observed brightness 
distribution is described by the expression: 



Jo 



2^lF 



exp 



Qgs - F{^)f 
2*2 



h(t- i p (V0) • 



• (t - t p (ip)^ exp 
where t p (if) = t p tjj s = F(tp)j . 



-+ \ 2 



2T S 



(t-tp($) 



(21) 
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Normalized quantities v = tp/^i, v s = * s /^, vo = ^o/^i, F(v) = F(* z i7)/* z , 
and normalized time intervals r = t/T s and t p {u) = t p {^>iP, ^>iF(V))/T s are more 
convenient to analyze gravitational focusing. In the normalized coordinates, the 
source brightness distribution (21) is reduced to a form: 



(22) 

We can get the light curve of the observed image J p {t) having integrated 
Eq. (22) over all possible angles of observations: 

/oo 
I p (v,r)dv. (23) 
-oo 

The full energy of the observed splash of brightness (impulse hereafter) can be 
determined as 



Jp 



/oo 
dr)J p (r). (24) 
-oo 



5. Calculation of light curves 

For the GL systems observed in optics, uq « 1 holds as a rule, excluding rather 
extended radio sources, where may be of the order of unity or even larger. 
It is known from the theory of GL, that the focusing effect of the gravitational 
field is negligible for the sources with vq > 1, therefore, this case will not be 
analyzed further. According to (21), for v$ « 1, radiation gets into the point 
of observation only from the small regions of space, situated near the point 
V = z/j, where the index of power in the exponent becomes zero, i.e. when 
u 8 = F(ui). The roots V = Ui of this equation, - lens equation (3), - determine 
the coordinates of brightness peaks in the images, formed by the GL. Supposing, 
that iV separate images of the source are formed by the lens, - i = 1, 2, . . . , N, - 
let us introduce small deviations Ai? from maxima Pi, - V = z?j + Av, - and make 
the following simplifications, (Minakov et al. 2001): 

V s - F{Vi + Az7) w (AuV)F(ui), 

t - T p (Pi + Az7) w T 'i - (Az7V)t p (z7;), (25) 

where t\ = r — T p {pi). With simplifications (25) taken into account, the dis- 
tribution of brightness (13) near the i-th image maximum may be presented 
as 



I pi (Au, T'i 



Jo 



1-Kvl 



h(r'i - AiAv x - BiAvy) ■ (r'j - MAv x - BiAv y )x 



x exp 



x exp 



--(r'i- AiAv x - BiAv y ) 



lf~ ~ ~ ~ \: 

^2 \ {CiAv x + DiAv y )e x + {DiAv x + EiAv y )e y } 



(26) 
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where 

A - 9t p. r - 9t p. r - dFx - n - dFy - p - dFx - dFy e>7\ 
az^ af^ ou x ou y ou y ov x 

All the derivatives in (27) are calculated in the point V = z?j. Having in mind the 
introduced simplifications, after simple calculations, one can get the following 
asymptotic approximation for the light curve of the i-th image: 



/oo 
d(Av)I pi {AiJ, r'i) 
-oo 



' l + Ri ' eXP \ 2 ' Ti 



1 + v 7 ^ • . , = ■ exp ' l — . 1 + $ ' 



(28) 



y/2(l +Ri) *\2(l + Ri)J \ \y/2(l + Ri)J, 

The following designations are introduced here: 

R i = "irUr! R i = (C? + Df)Bl + (A 2 + El)~A\ - iMB&XCi + Ei); (29) 

q i = \C i E i -D 2 i \-\ (30) 

and $>(x) is the probability integral. The introduced quantity qi can be easily 
shown to represent the amplification factor for the i-th image of a stationary 
source. 

As it follows from simple physical considerations, amplification of the ini- 
tial impulse Jo in gravitational focusing originates only from the spatial redis- 
tribution of brightness, determined by the lens amplification factor q. Thus, a 
condition J p i = qiJo must hold. Indeed, having integrated (28) over the whole 
region of variation of r'i, (— oo < r'j < oo), we get the full energy of the focused 
impulse, arrived to the point of observations from the i-th image: 

/oo 
J P i{T'i)dT'i = qiJ . (31) 
-oo 

For the values Ri » 1, we get the following approximate expression for 
Jpi(T~'i), having used the asymptotic expression of <J>(x): 

Jpiir'i) = qi J • r'i • exp (^^fj . ^ i > 0- (32) 

The light curve of the observed i-th image is seen to be formed from the 
initial undisturbed light curve 

4°\t') = J .TVex P (-^V t\>0 (33) 



by multiplying the latter by the amplification factor It should be noted 
however, that such a simple effect of a GL is observed only for i?j >> 1. At small 
values of Ri, the lens will distort the initial light curve of the source. In addition, 



12 Minakov k Vakulik 



the amount of distortions will be different for different images because of the 
existing discrepancy in the values of Ri. Relationships Q{t') = \/2tt J p (rf)/ Jog 
for different R are presented in Fig. 7. 

We investigated focusing of a quasar radiation by a regular component of 
the Q2237+0305 system as an example. For this system, F(p) = Fm{v) and 
is described by Eq. (4), while the signal time delay r p (z/), disregarding an 
insignificant additive term, is described by expression: 



t p {v) = t + t< 



2Z 



(1 + a 



2 V 2 - 2 + 



1 



+ 2a 1 - 



1 



(-1 ~ -I) 



a 



(34) 



where T g = 2r g /cT s , r g ~ 3 • 10 9 -r- 3 • 10 10 /cm is a gravitational radius of the 
galaxy nucleus, and Z s /Z « 10 is the relative distance, (Kayser, Refsdal & 
Stabell 1986). It is easy to show, that there are four points in Eq.(34), where 
T p (u) is minimal. For a = 0.16 and Z s /Z = 10, coordinates of these points are, 
respectively: 



^1,3 



Z a (l + ct)-Z 



0.97, 991 = 0, ip 3 = tt; 



V2A 



Z s (l -o)-Z 



1.16, 



<^2 



7T 

2' 



V? 4 



3vr 
T' 



(35) 



Table 1. Values of normalized angles and parameters for 
Q2237+0305 model 



Image 


v x 


Vy 


Vx{y 8 = 0) 


Vy(V=Q) 


1i 


7* 


ft 


A 


0.325 


-1.031 





-1.08 


2.261 


0.251 


0.446 


B 


0.279 


1.041 





1.08 


2.244 


0.258 


0.443 


C 


-0.88 


-0.044 


-0.94 





0.969 


5.21 


0.774 


D 


0.994 


-0.099 


0.94 





3.395 


0.015 


0.526 


Nucleus 





















Presense of minima in t p {v) indicates, that temporal changes in the source 
will be observed in the lensed images in the directions of points (35) at first, and 
then in the A, B, C and D images. 

For the obtained t p (0), (34), a dimensionless parameter Hi in (29) can be 
represented as Hi = Tgj(yi), where = ji = Hi/i~g does not depend on 

T g and determines only by the value of a, relative distance Z s /Z, and by the 
coordinates of the four images. For the given coordinates of the images, the cor- 
responding values of parameters 7, and ft = [r p (^) — to]/t 9 , (i = A,B,C, D) 
were calculated. The results are presented in Table 1. The coordinates of the 
observed lensed images from (Crane et al. 1991) are also presented here, as well 
as the coordinates of the same images for the central source {y s = 0). 
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Further calculations were made for the angular dimension of the source of 
v k, 10~ 3 , and for two values of r g = 0.1, r g = 600. In Fig. 7, the individual light 
curves J p i{r' \) (left), and a total one J p s = Jpii^'pi) f° r T g = 0-1) (right) are 
presented. 

It is evident, that there is no possibility to detect separate images in the 
total light curve at small values of r g . The components become distinguishable 
in the total light curve for sufficiently large values of T g . For example, a total 
light curve for T g = 600 is shown in Fig. 8. 

The discussed light curves referred only to a variable component of the 
source radiation. In reality, the radiation can be represented as a sum of a 
stationary and variable constituent. In this case, with (31) taken into account, 
the initial light curve would be given by the following expression: 



J$(T') = J s + JoT'ex P {--T>i), (36) 

where J s = const is a stationary component of the source radiation. The fol- 
lowing total brightness as a function of time will be recorded in the point of 
observation: 



J P T.{r') = qJ s + Jo l^'i ex P i y ~\ T> ^) = Js q + Y^- / qiT>i 6XP (~^ T '* 2 



2 

where q = ^ Qi ~ 8.9 is a total amplification factor of all the four images of 
Q2237+0305. 



6. Discussion and conclusions 

Some conclusions can be made from the analysis of gravitational focusing of the 
extended source, that varies in time. 

If a gravitational lens forms several images of the source, the light curve 
of an individual image J p (r') will depend on the value of parameter R, which 
accounts simultaneously for geometry of the task (distances from the observer 
to the source and to the lens, Z p and Z s respectively), parameters of the lens 

(a total mass and its distribution), and the source parameters (\I/o> ^s,Ts)- For 
the model of GL Q2237+0305 under consideration, - a lens with a quadrupole, 
- the value of R was estimated to be R ~ {yo T g) ~ 2 ■ 

If R » 1, {voT g « 1 - a "small" source, "weak" lens or "long" impulse), 
the light curve J p (t') does not depend on R any more, and differs from the 

undisturbed light curve J^\t') by a multiplier q only. 

For small values of R, when R « 1, {vQT g > 1 "extended" source, 
"strong" lens, or "short" impulse), a smoothing action of the lens should be taken 
into account. Amounts of deformations will be different for different images, 
because of the existing difference in R. 

The analysis of focusing action of a GL with an arbitrary mass distribution 
allows to argue, that the GL acts as a filter, transmitting low frequencies (slow 
brightness variations) with no distortion, and "rejecting" the higher ones, - rapid 
variations. 
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Specific durations of the source brightness variations, T s , for which a smooth- 
ing action of a GL should be taken into accout, can be determined as T s < r g vo/c, 
(vqTq < 1). For example, for a regular model of Q2237+0305, supposing the mass 
of the galaxy nucleus M ~ 1O 1O M0, (r g ~ 3 • 10 10 /cm), and choosing the charac- 
teristic angular size of the source to be equal to vq ~ 10~ 4 , we get the estimate 
for T s < 10s. Observations of the Einstein Cross in the IR spectral range, (Agol, 
Jones k Blaes 2000), allow to estimate vq as 10" 5 < v o < 10 -1 . Having adopted, 
e.g., z^o ~ 10~ 2 , we get T s < 10 3 sec. 
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Figure captions 

Fig 1. Mutual locations of a source, S, the observer P, a macrolens GL, and a 
microlens m. 

Fig 2. Critical curves (1, 2, 3) and caustics (l',2',3') for a microlens location a) 
outside, b) inside and c) exactly at the critical curve of the macrolens (3). The 
curves (1, 1') are for the microlens displacement along the X axis, and (2, 2') - 
along the Y axis; 3'is a macrocaustic. 
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Fig. 3. Deformations of a macrolens critical curve Av(ip) = v cr {ip) — v££(<p) for 
different distances v m = (v m ,(Pm = 7r) between a microlens and the macrolens 
critical curve. 

Fig 4. Critical curves (left) and caustics (right) of the discussed macrolens 
model of Q2237+0305, with 300 microlenses, distributed homogeneously within 
the strip of 2a width along the macrolens critical curve v^J. . 

Fig 5. Difference of the critical curves (2) from those calculated in the linear 
approximation (1), for V = 10 -9 , and for the microlens locations a)outside, 
b)inside, and c)exactly at the critical curve of the macrolens (3). The microlens 
coordinates are (0, 0). 

Fig 6. Behavior of Q(t') = \/2ttJ p (t') / J^q at different values of R. 

Fig 7. Individual (left) and summarized (right) light curves for the values of 
parameter r g = 0.1. 

Fig 8. A total light curve for the value of r g = 600. 
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